Magnetizing a complex plasma without a magnetic field 
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We propose and demonstrate a concept that mimics the magnetization of the heavy dust particles 
in a complex plasma while leaving the properties of the light species practically unaffected. It makes 
use of the frictional coupling between a complex plasma and the neutral gas, which allows to transfer 
angular momentum from a rotating gas column to a well-controlled rotation of the dust cloud. This 
induces a Coriolis force that acts exactly as the Lorentz force in a magnetic field. Experimental 
normal mode measurements for a small dust cluster with four particles show excellent agreement 
with theoretical predictions for a magnetized plasma. 

PACS numbers: 52.27.Gr, 52.27.Lw, 52.25.Xz, 45.50.Jf 
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FIG. 1. (Color online) Estimates of u c /u)o in various strongly 
coupled plasmas: white dwarf stars [7] (p m > 10 9 kgm -3 , 
B < 10 5 T), neutron stars [14] (p m ~ 10 7 . . . 10 14 kg m" 3 , B < 
10 11 T), ions in Penning traps [Hi [TJ] (n ~ 10 7 . . . 10 8 cm" 3 , 
fields of several Tesla) , and complex plasmas [8] 1101 111] . The 
method presented in this paper (rotating complex plasmas) 
covers a wide range of magnetizations at room temperature. 

Strongly coupled Coulomb systems have very unusual 
properties including spontaneous spatial ordering and the 
formation of liquids or even crystals. Coulomb crystals 
and liquids, originally predicted by Wigner pQ, were even- 
tually observed on the surface of helium droplets , in 
ion traps [3], and in complex plasmas 4, 5 . They are also 
believed to occur in semiconductor quantum dots and 
quantum wells [S] , as well as in white dwarf and neutron 
stars [7J. Of particular current interest is their behavior 
in a magnetic field, where strongly modified oscillation 
spectra [8HTT] or anomalous diffusion properties [12] have 
been predicted, and even applications to verify enhanced 
nuclear reaction rates have been demonstrated |13j . Es- 
pecially in neutron stars |14j . giant magnetic fields are 
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present that considerably modify the properties of the 
liquid and crystal states in the outer layers and alter the 
whole evolution of the star. 

In a strongly coupled plasma (SCP), the Coulomb in- 
teraction energy of two particles, Q 2 / (Aneoa) [charge Q, 
typical inter-particle distance a], is much larger than 
their thermal energy, ksT . In a magnetized SCP the 
particles are, in addition to the electrostatic interac- 
tions, subject to the Lorentz force. The degree of mag- 
netization can be measured by comparing the relevant 
time scales associated with these two forces [IB]: While 
the Coulomb interactions lead to a characteristic vibra- 
tion frequency [S] ojq = y / Q 2 /(47re ma 3 ), the cyclotron 
frequency uj c = QB/m is the relevant parameter for 
the Lorentz force. For a given magnetic field B and 
mass density of the plasma p m = urn (particle mass 
to, number density n), the ratio of the two becomes 
u; c /cl>o = By^Seo/ p nl . Estimates for various strongly cou- 
pled astrophysical [7j [14] and laboratory plasmas (T3l [15] 
are presented in Fig. [I] 

Complex (dusty) plasmas |17| have, in recent years, 
become a prototypical system to study strong correlation 
effects in unmagnetized Coulomb systems. They con- 
tain highly charged, micrometer-sized particles embed- 
ded in a partially ionized electron-ion plasma. Complex 
plasmas are found in numerous space environments in- 
cluding interstellar clouds, cometary tails, or planetary 
rings [18] . In laboratory experiments, the large particle 
size and mass make it possible to follow individual parti- 
cle trajectories with unprecedented spatial and temporal 
resolution [T§H2"T] providing valuable insight into strong 
coupling phenomena. However, until now, it has not been 
possible to extend this analysis to magnetized Coulomb 
systems. The large particle mass (via large p m ) limits 
uj c /u>o to values below 0.1... 0.5, even if superconduct- 
ing magnets and particles with sub-micron diameter are 
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used [HI HO] j see Fig. [I] Also, inevitably, such fields rad- 
ically alter the plasma conditions and may lead to fila- 
mentation :22J. 

Here, we propose a completely different approach to 
achieve very strong magnetization effects that overcomes 
all these limitations. It can be used with standard 
(micron-sized) particles, does not significantly alter the 
plasma properties, and does not require expensive mag- 
nets at all. The basic idea to effectively "magnetize" the 
complex plasma is to impose the global motion of a ro- 
tating neutral gas column on the dust particles, which 
can be achieved, e.g., by a rotating electrode (23] ■ I n 
the frame co-rotating with the neutral gas, the Coriolis 
force takes over the role of the Lorentz force and mimics 
the effect of a magnetic field, while the centrifugal force 
renormalizes the confinement. In all other aspects, the 
equations of motion are equivalent to those in a gas at 
rest. Related methods, where rotation is induced by dif- 
ferent means, have been applied in the context of cold 
quantum gases [2"1H2"7] . 

In the following, we first present the theory underlying 
our concept and work out the scaling relations for the 
plasma parameters in a two-dimensional harmonic trap. 
We then perform a proof-of-principle experiment for the 
normal modes of a small cluster and find excellent agree- 
ment with the theoretical predictions. 

In our model, the neutral gas is assumed to rotate as a 
rigid body with a velocity profile u(r) = (il e z ) x r. Here, 
e z denotes the unit vector in the z-direction and fl the 
rotation frequency. Particle confinement is provided by a 
harmonic trapping potential V(p, z) = ? ( w j_P 2 + u 2 ? 2 ) > 
where p = \J ' x 2 + y 2 . The equations of motion for the 
dust particles then have the form (i = 1, . . . , N) 

N 

mr, = - ViV( Pi , Zi ) + Fjf ~ vm [n - u(n)] + f t . 

Here, F-j* = —\7i4>(pij, z^) denotes the interaction force 
[zij = Zi — Zj, pij — + yfj], 4>{p,z) the associated 

potential [28] . v the dust-neutral friction coefficient, and 
fi(t) the stochastic force. The latter has zero mean and 
the correlation (/f (*)//(*')) = 2mvk B T n 8 lJ 5 a ' 3 5{t - t'), 
where i,j € {1, . . . , N} and a, /3 € {x, y, z}. The neutral 
gas temperature is denoted by T n . 

In a frame rotating around the z-axis at the angular 
velocity ft of the neutral gas, the equations of motion 
take a form that illustrates the influence of the gas flow. 
Denoting the coordinates in the rotating frame with an 
overbar, the relation between the laboratory frame and 
the rotating frame reads r(t) = R(t)f(t), where R(t) 
is the matrix for rotations around the z-axis (see, e.g., 
Ref. [1H])- Replacing the coordinates leads to 

N 

mh = -ViV{pi, Zi) + Fif + FcoAh) ~ vrnfi + f l: 



where the effective confinement potential is now given by 

HM) = f [K-o 2 )p 2 + ^ 2 z- 2 ]. 

The first observation is that, in the rotating frame, 
the effective confinement frequency Q± = \J iq\ — Q 2 in 
the direction perpendicular to the rotation axis is being 
reduced as a consequence of the centrifugal force. Second, 
the rotation induces the Coriolis force 

-Fbor(^) = mf X (20 e z ), 

which has the same form as the Lorentz force for a homo- 
geneous magnetic field B cS — (2mQ/Q)e z . The doubled 
rotation frequency can be identified with the cyclotron 
frequency, uj c — 251. The neutral gas flow velocity u(r) 
no longer appears explicitly because the coordinate frame 
rotates at the same angular velocity as the gas. While 
the form of the inter-particle forces remains unaffected, 
-F™* = — V i<t>{p~ij , z^) , the components of the random 
force occur in a mixed form, fi(t) = R T (t)fi(t). How- 
ever, its statistical properties are the same as in the lab- 
oratory frame, i.e., a Gaussian white noise is retained, 
(fi(t)) = 0, and (f?(t)fP(t')} = 2m V k B T a S lJ S a ^S(t~t / ). 

Let us now estimate the effective "magnetic fields" 
|-Beff| = 2mfl/Q that can be reached. With rotation 
frequencies of 51 « 10 Hz, a charge of |Q| « 10 4 e, and 
a mass of m w 10 -12 kg, one easily generates effective 
magnetic fields exceeding 10 4 T, which is far beyond the 
capabilities of superconducting magnets. At very high 
rotation speeds, the centrifugal force weakens the hor- 
izontal confinement j3IJ|, which may give rise to a de- 
formed plasma shape and a lower dust density. However, 
especially in two-dimensional systems [31] . the shape of 
the dust monolayer will remain unaffected, and a lower 
density could even be beneficial as it decreases loq, thus 
allowing to reach even larger magnetizations 2CI/luq. 

To illustrate these effects in more detail, we con- 
sider the limit u z luj_. This situation is typi- 
cally encountered in dusty plasma experiments in radio- 
frequency discharges, where the particles are located 
in the sheath region above the lower electrode. Their 
interaction is well described by the Yukawa potential 
4>(r) = (Q 2 /4ireo r) exp(— r/A) with the Debye screening 
length A. As a typical inter-particle distance we choose 
a(Q) = [Q 2 /(47reo mu})] 1 ' 3 , which is equivalent to (up 
to a numerical factor) the ground state distance of two 
trapped particles with Coulomb interaction. The effec- 
tive trap frequency u>j_ then coincides with the previous 
definition of ujo = \J Q 2 /(4:neo ma 3 ) and provides a nat- 
ural frequency unit for particle vibrations. Both the dis- 
tance and time scales, a(il) and wj 1 (0), increase upon 
faster rotation. Since the electron and ion time scales are 
orders of magnitude larger than the rotation frequencies, 
the charge Q and the Debye length A (i.e., the plasma 
properties) remain constant to high accuracy. 

The problem is then governed by the following di- 
mensionless parameters: the coupling parameter T(il) = 
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FIG. 2. (Color online) Left: Variation of the effective mag- 
netization, the coupling parameter, screening parameter and 
damping constant with the rotation frequency in a finite two- 
dimensional system. Right: Scaling of F and function of 
the effective magnetization 2fi /ui± for typical complex plasma 
conditions. 

Q 2 /(47re a(i7) fc B T n ), the screening parameter = 
a(f2)/A, the damping rate 7(f2) = v/uj±(il), and the ef- 
fective magnetization 

Wx(fi) (1 -fi 2 /w2)l/2- 1 j 

It is obvious that, already for a relatively slow rotation 
with fl = oj±/2, the plasma is strongly "magnetized" 
(2fi/cDj_ = 1-15), see the left panel of Fig. [2] The asso- 
ciated changes (compared to the non-rotating system) of 
the coupling and screening parameter are small. These 
parameters scale as r(O)/F(0) = (1 - Q, 2 /cj^) 1 / 3 and 
k(O)/k(0) = (1 — fi 2 /^) -1 / 3 , respectively, and start to 
change substantially when the centrifugal force notice- 
ably increases the inter-particle distance (Q/uj± > 0.7), 
see the left panel of Fig. [2j In this regime, the decrease 
of Q±_ is largely responsible for the dramatic increase 
of the magnetization parameter. We also show the di- 
mensionless damping rate, which scales as 7(f2)/7(0) = 
(1 — n 2 /lj 2 )^ 1 / 2 . This means that, in experiments, the 
gas pressure should be sufficiently low allowing for a small 
neutral gas friction coefficient v before start of the rota- 
tion. The right panel of Fig.[2]shows T and function 
of the effective magnetization in a parameter regime that 
should be easily accessible in dusty plasma experiments, 
see below [52"] , 

In the following, we present a proof-of-principle exper- 
iment to verify the efficiency of the proposed concept. A 
sketch of the experimental setup is shown in Fig. [3| The 
experiments were performed in a 13.56 MHz capacitively 
coupled radio-frequency discharge at a gas pressure of 
p = 0.4 Pa (Argon). Spherical particles with a diameter 
of d = 21.8 /im and a mass of m = 6.46 • 10~ 12 kg are in- 
jected into the plasma, where they form two-dimensional 
clusters. The upper electrode can be set into rotation 
with frequencies up to 30 Hz, which causes a vertically 
sheared rotational motion of the neutral gas column [53] . 




FIG. 3. (Color online) Experimental setup: Dust particles 
are confined in a layer approximately 1 cm above the driven 
electrode. The electric field within the plasma sheath com- 
pensates gravity and provides vertical confinement, while hor- 
izontal confinement is realized by a cylindrical cavity of 20 mm 
in diameter and 2 mm in depth. The particles are illuminated 
from the side by a laser fan and can be observed from the bot- 
tom through a transparent electrode coated with indium tin 
oxide (ITO). The top electrode is rotatable. 



We concentrate on the dynamics of a small ensemble 
of N = 4 particles and analyze their normal modes. The 
normal modes of small 2D clusters have already been 
studied experimentally 3.Tj. but for magnetized dusty 
plasmas only theoretical predictions exist [34] . We first 
consider the center-of-mass (sloshing) mode. Since the ef- 
fective confinement in the rotating frame is harmonic, the 
center-of-mass coordinate f cm (t) is independent of the in- 
teraction (Kohn theorem) and obeys the same equation 
of motion as a single particle [35 . In the absence of rota- 
tion, the two center-of-mass modes are degenerate with 
uj cm = uj_\_. For f2 > 0, however, this degeneracy is lifted, 

and the frequencies read uf m — \J (lo c /2) 2 + lj 2 _±uj c /2 = 

luj_ ± n. Here, uj± = — fl 2 is the effective trap fre- 

quency in the rotating frame. The experimental results 
obtained from the spectrum of f cm (t) are depicted in 
Fig. [4] and show remarkable agreement with the theoret- 
ical prediction. 

It is crucial for our scheme to further verify the accu- 
racy of the remaining modes, which are sensitive to both 
the magnetic field and inter-particle correlations. By lin- 
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earizing the equations of motion in the rotating frame, 
we determine the eigenfrequencies u> and eigenvectors i>, 

0, 
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where a, /3 £ {x, y} and i, j £ {1, . . . , N}. Further, H"- p 
denotes the Hessian of the total potential energy in the 
rotating frame, Sij (S a/3 ) the Kronecker delta, and e a @ z 
the Levi-Civita symbol. The cluster configuration for 
N = 4 is a square with particles located a distance R 
from the trap center, which is calculated from Eq. (6) 
in Ref. Even though a finite dust-neutral friction 

parameter is essential to put the plasma into rotation, it 
is sufficiently low to be negligible for the calculation of 
the eigenfrequencies (v/uj± s» 1/50). To determine the 
mode frequencies experimentally, we calculated the pro- 
jection P(t) = J2iLi ' v i 01 the particle trajectories 
f i (t) on the eigenvectors Vi of the unmagnetized system, 
see Fig. |4j The breathing mode corresponds to a radial 
expansion and contraction of the cluster. The spectrum 
of P(t) shows a peak at the associated mode frequency, 
which can be tracked as the rotation frequency is varied. 

The measurements are compared with the theoretical 
results in Fig. [4] As for the case of the sloshing modes, we 
observe excellent agreement with the normal modes of a 
magnetized plasma. Effective magnetizations 2f2/a/j_ > 3 
allow us to clearly verify the predicted splitting of the 
normal modes into the upper and lower branch [34] . 

To summarize, we have presented a simple approach to 
"magnetize" a complex plasma. The idea is based on the 
correspondence between charged particles in a magnetic 
field and particles in a rotating gas column. The pos- 
sibility to put the plasma into rotation takes advantage 
of the dissipative nature of complex plasmas in which 
the neutral gas acts as a highly effective transmission 
agent of angular momentum. We demonstrated that, 
with very modest rotation frequencies applied to only 
one electrode, strongly correlated particles in a rotating 
dusty plasma reproduce, to high accuracy, the normal 
modes of a magnetized system. Evidently, our concept 
can also be realized by other means and opens new unique 
possibilities for highly accurate studies of strongly corre- 
lated and strongly magnetized plasmas. The advantage 
compared to cryogenic ions is the broad range of accessi- 
ble plasma parameters that can be varied independently 
(coupling strength, screening, dissipation) and the avail- 
ability of single-particle resolution at room temperature 
conditions. It should be possible to create SCP states 
with extreme magnetizations potentially even compara- 
ble to those exotic ones in the outer layers of neutron 
stars. While our technique does not require any (super- 
conducting) magnet at all, use of the latter in combina- 
tion with plasma rotation allows to create novel types of 
plasmas. In such plasmas, there would be effectively two 
"magnetic fields" that can be controlled independently — 
one of which (the rotation) affects only the heavy parti- 
cles whereas the second ( "real" ) field influences the clec- 




FIG. 4. (Color online) Comparison of experiment and theory: 
Frequencies of the seven non-trivial normal modes of N = 4 
dust particles as a function of the scaled rotation frequency. 
Symbols denote experimental results (uj±/2ir — 2.52 s _1 ), 
lines are the theoretical eigenfrequencies with k(0) = 0.7. 
The eigenvectors in the unmagnetized limit 36 (fi — 0) are 
sketched for four particular modes. Note that the two modes 
(iii) are degenerate in this case. The breathing mode and the 
two center-of-mass modes are indicated by "br" and "cm*", 
respectively. The effective magnetization 2Q/oj± is shown on 
the upper axis [Eq. Q]. 



trons and ions, allowing for an effective control of the 
interaction between the dust particles. The combination 
of a rotating flow with a magnetic field also provides a 
promising avenue to access magnetorotational instabili- 
ties on the particle scale, which are relevant for the un- 
derstanding of accretion disks [37]. 

This work was supported by the DFG via SFB TR6 
and SFB TR24. 
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